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Abstract. The amount of late decaying massive particles (e.g., gravitinos, moduli) produced 
in the evaporation of primordial black holes (PBHs) of mass Mbh ^ 10^ g is calculated. 
Limits imposed by big-bang nucleosynthesis on the abundance of these particles are used to 
constrain the initial PBH mass fraction f3 (ratio of PBH energy density to critical energy 
density at formation), as: /3 < 5 x 10-^^{x^/Q10-^)-\Mbh/}0^ g)~^^'^(Ki,/lO-^^); is the 
fraction of PBH luminosity going into gravitinos or moduli, is the upper bound imposed 
by nucleosynthesis on the number density to entropy density ratio of gravitinos or moduli. 
This notably implies that such PBHs should never come to dominate the cosmic energy 
density. 
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1. Introduction - The spectrum of locally supersymmetric theories generically contain 
fields whose interactions are gravitational, and whose mass ~ (9(100 GeV). The Polonyi 
and gravitino fields of supergravity theories, or the moduli of string theories, are such ex- 
amples. This leads to well-known cosmological difficulties: quite notably, such particles 
(hereafter generically noted and termed moduli) decay on a timescale ~ Mpi/m^ ~ 
10^ s (m^/100 GeV)~^, i.e., after big-bang nucleosynthesis (BBN), and the decay products 
may drastically alter the light elements abundances 0]. The success of BBN predictions 
provides in turn a stringent upper limit on the number density to entropy density ratio (Y"^) 
of these moduh, generically < 10~^^ 0] (see Sec. 3). 

It is argued in this letter that these same constraints can be translated into stringent con- 
straints on the abundance of primordial black holes (PBHs) with mass Mbh ^ 10^ g. In effect, 
moduli are expected to be part of the Hawking radiation of an evaporating black hole as soon 
as the temperature of the black hole exceeds (roughly speaking) the rest-mass m^; and indeed, 
the Hawking temperature of a PBH reads Tbh = m^jMBH ~ 10^ GeV {Mbr/IO^ g)'^ §. 

Primordial black holes are liable to form in the early Universe at various epochs, e.g., 
when a density fluctuation re-enters the horizon with an overdensity of order unity [^], or 
when the speed of sound vanishes (as may occur in phase transitions). As a consequence, 
constraints on the abundance of PBHs can be translated into constraints on the structure 
of the very early Universe Q. Until recently, the only existing constraint on PBHs of mass 
^BH ^ 10^ g relied on the assumption that via evaporation, PBHs leave behind stable Planck 
mass relics . However, recent work from the perspective of string theories seems to indicate 
that this is not the case in particular that evaporation proceeds fully. Nevertheless, 
Green |^ has pointed out recently that such PBHs would also produce supersymmetric 
particles, and consequently, cosmological constraints on the lightest supersymmetric particle 
(LSP) density could be turned into constraints on the initial PBH mass fraction f3 (defined 
as the ratio of PBH energy density to critical energy density at formation). This constraint 
relies on the assumption that the LSP is stable, i.e. i?— parity is a valid symmetry; and, as 
attractive as i?— parity is, it is not of a vital necessity altogether. The constraint related 
to the production of gravitinos or moduli, to be derived below, is thus complementary to 
this -R— parity constraint, and it also turns out to be more stringent. Hereafter, units are 
h = k-B = c = 1, and mpi = Mpi/(87r)^/^ ~ 2.4 x 10^*^ GeV is the reduced Planck mass. 

2. Moduli production - Although one is generally interested in itself, and not in its 
momentum dependence, it will prove necessary in a first approach to keep track of dY^/dk 
(where k is the momentum) integrated over the black hole lifetime. In effect, during their 
evaporation, PBHs produce moduli over a whole spectrum of momenta, with high Lorentz 
factors, and the existing constraints on depend strongly on the (cosmic) time at which 
moduli decay (r^ is the decay timescale in the modulus rest frame), hence on whether they 
are relativistic or not. 

More quantitatively, the mass and temperature of a PBH evolve with time t during 
evaporation as: M{t) = Mbh [I - (t - ti)/TBnf^ and T{t) = Tbh [I - (t - ti)/rBH]"'^' l- 
Here, ti denotes the time of formation, ti <^ tbh, with tbh the PBH lifetime: tbh — 
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0.14 s(Mbh/10^ g)^ Q Toward the end of the evaporation process, the temperature increases 
without apparent bound, although the standard analysis breaks down at T ~ mpi (see 
Ref. for a discussion of the end point of evaporation). Once the black hole temperature 
T 3> m^, moduli can be considered as massless. Then the number of moduli emitted per 
PBH, with momentum k between k and k + dk, and per unit of time, is, for a Schwarzschild 
black hole ||: q4k,t) = {2ny^T^{M{t), k)/ [exp{k/T{t)) - {-if']. The absorption coeffi- 
cient is a non-trivial function of M, k and s which has to be calculated numerically [|T0| , 



and s is the spin of (j). As announced any PBH will thus produce moduli at some point, 
and, moreover, these moduli will be produced over a whole range in momentum. To give 
an example of the sensitivity of the constraints on on the time of decay: if decays into 
photons, pair creation on the cosmic background (of temperature T^) suppresses cascade 
photons whose energy E > m1/22T^] since ~ lMeV(t/ls)~^/2^ at early times < lOS, 
the cut-off lies below the threshold of deuterium photo-dissociation (~ 2 MeV), and the con- 
straints on are evaded, while at later times, the cut-off is pushed above this threshold, and 
photo-dissociation becomes highly effective. Finally, since a modulus carrying momentum k 
at cosmic time will decay at time t ~ T^m.a,'x.[{k / m^Yf^ , 1], it is necessary to follow dY^/dk 
as a function of time. As an aside, this will permit the calculation of produced by PBHs 
such that Tbh < "^0- 

This calculation is carried out below in the following limits. As a first approximation, it is 
sufficient to assume that all particles are emitted at the same average energy, parametrized 
as aT{t); a is a constant which depends on s, with a ~ 2.8 for s = 0, a ~ 4 for s = 1/2, 
and a ~ 7 — 8 for s = 3/2 [jll|.This approximation suffices as the energy at peak flux 
corresponds to the average energy to within ~ 10% []TT|, and since the injection spectrum 



cuts-off exponentially for k > aT, and as a power-law for k < aT. The initial mass fraction of 
PBHs is approximated to a delta function centered on Mbh- Although recent considerations 



tend to indicate otherwise |]T2[, this remains a standard and simple approximation; moreover, 
the extension of the results to a more evolved mass fraction is easy to carry out. Finally, 
it is also assumed that the Universe is radiation dominated all throughout the evaporation 
process, which implicitly implies that black holes never dominate the energy density. This 
latter assumption will be justified in Section 3. 

Then the distribution f^{k,t) = s^^dn^/dk = dY^p/dk, where s denotes the radiation 
entropy density, at times tbh < t < r^f, reads: 

■'^H dk' 



U{k,t)=Y^nl^%4k',t')—dt'. 



t The lifetime of a black hole depends on the number of degrees of freedom gs in each spin s in the 
radiation [[l), i.e. tbh = 6.2s/(MBH)~^(A^BH/10^g)^ with /(Mbh) ^ 0.267go + 0.147.gi/2 + O.O651 + 
0.02(73/2 + 0.007(72- Here the particle content of the minimal supersymmetric standard model (MSSM) with 
unbroken supersymmetry has been used, go = 98, 51/2 = 122, gi = 24, 513/2 = 2, 52 = 2. 



•tXhe value of a for s = 3/2 is based on extrapolation of the results of Ref. |11 for other spins, while the 
fraction of luminosity emitted in spin s = 3/2 (noted in the following) is given in Ref. |^]. It docs not 
seem that a detailed study of Hawking radiation of gravitinos has ever been performed. Here it is assumed 
that the helicity states ±1/2 and ±3/2 of the gravitino are produced with values of a and as quoted for 
generic spin s = 1/2 and s = 3/2 respectively. 
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In this expression, q^{k',t') is the injection spectrum per black hole as above, Ybh = 
^bh/s, where tzbh represents the PBH number density, and k' = ka{t)/a{t'), where a is 
the scale factor. The factor dk'/dk results from redshifting of k' at injection time t' down 
to k at time t. Equation can be derived as the solution of the transport equation: 
dtfif, = Hdk{kf^) + Y-Bnqik.t), where H is the Hubble scale at time t, and the first term 
on the r.h.s accounts for redshift losses. This equation and its solution Eq. (|l]) are valid for 
t <^ T^, when the decay of (p particles can be neglected. It should be recalled that in the 
range of masses m<^ and Mbh considered, indeed tbh ^ t,^- Equation (0) also neglects the 
entropy injected in the plasma by PBH evaporation, which remains a good approximation 
as long as PBHs carry only a small fraction of the total energy density at all times. 

For mono-energetic injection at k' = aT{t'): 

6[k' - aT{t')]. (2) 

Here denotes the fraction of PBH luminosity \dM/dt'\ carried away by moduli; for the 
MSSM content, x,^ ~ 6 x 10~^ for s = with one degree of freedom (e.g., a modulus field), 

^ 6 X 10^'^ for s = 1/2 with two degrees of freedom (e.g., helicity ±1/2 states of the 
gravitino), and ^ 9 x 10~^ for s = 3/2 with 2 degrees of freedom (e.g., helicity ±3/2 states 
of the gravitino) [ill (see also previous footnote). The S distribution can be rewritten as a 
function of t, using the identity: S[fit)] = \df /dt\'^6{t — tg), where tg is such that f{ts) = 
(here tg is uniquely and implicitly defined in terms of k, k'). Equation (|l|) can be integrated 
in the limits k <^ ko and k ^ ko, where kg = aTBuit / tbu)'^^"^ is the momentum at time t 
of a particle injected at time tbh with momentum oTbh- In particular, modes with k <^ ko 
were injected with energy ~ aT^n at time t' ~ t{k/aTB}iY <^ tbh, while modes with k^ k^ 
were produced in the final stages at time t' ~ tbh with momentum k' ~ aT(t') ^ qTbr- 
One obtains: 



aT{t') 



dM 



dt' 



U{k,t) ~ ^x^ ^^" i; ^Fbh ik<^ko), (3) 

3 (oTbh)^ "Tbh Tbh 
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and both expressions agree to within a factor 3/2 at k = ko- 

If initially Tbh < ^(j>, moduli are produced only in the final stages for k' ^ at 
injection. Hence the above spectrum should remain valid if a low-momentum cut-off kc ~ 
^4>{t/TB}i)~^^'^ > ^0 is introduced. The total number of (f) particles produced (hence the 
constraint on /5) is thus suppressed (weakened) by a factor ~ {itT'^/TbrY, after integration 
of f^{k,t) over k > k^ if Tbh < ^(t>, i-e., if Mbh > 10^ g (m^/lO TeV)"^. Since this mass 
range Mbh ^ 10^ g is moreover strongly constrained by the effects on BBN of quarks directly 



produced in the evaporation |]T3|, it will be ignored in the following. 

For PBHs such that Tbh ^ it is a very good approximation to consider that emitted 
moduli carry at time t a momentum ko, since kf^{k,t) = dY^/d\n{k) behaves as k"^ for 
/c -C /co, and as /c~^ for k^ k^. Moreover at time t = t^: 
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and therefore the (j) particles decay at rest (in the plasma rest frame), at time ~ 
Mpi/m^, in the range of masses considered, < lOTeV and Mbh ^ 10^ g. One then seeks 
the total number of moduli present at that time, which is given by: 

Mbh ^ 

This result can be obtained as a solution of the transport equation 
dfY^ = x^YBu\dM{t)/dt\ /aT{t), or by integrating f^{k,t) over k in Eqs. (|^), above; 
all three results agree to within a factor 3/2. Equation has a simple interpretation: 
within a factor 2 it corresponds to the instantaneous evaporation of black holes at time 
tbh, with total conversion of their mass Mbh in particles of energy oTbh, a fraction 
of which is moduli. This result can be rewritten in terms of more conventional parame- 
ters. The mass Mbh is taken to be a fraction 5 of the mass within the horizon at the time 
of formation t^: Mbh ~ AnSrripi/ Hi, where Hi denotes the Hubble scale at time tj, and 
6 ~ 0{1) is expected [|T^. Furthermore, instead of Ibh, one generally uses the mass fraction 



/3 = n^nMBu/pc defined at the time of PBH formation ti, with = SHfrripi the critical 
energy density at that time. Using s = (27r^/45)5'^,T^, with the cosmic background tem- 
perature, ~ 0.5g20(I^Hl^'^mp-^ , and (7200 = fi'^/200 {g^, number of degrees of freedom), one 



obtains: 




/?^3xlO-,2^/oV^/^ ^ Fbh, (7) 




and therefore: 

^ 2 X lO'S'/'gJ' 

which constitutes the main result of this section. If the Universe went through a matter 
dominated era between times ta and th, with ti < ta < h <ti tbk, then the r.h.s. of Eq. (|) 
must be multiplied by the factor (Hf,/ HaY^'^, where Ha,b is the Hubble scale at time tafi, and 
the constraint on /? Eq. (|) below is weakened consequently. 

3. Discussion - As mentioned previously, the most stringent constraints on result 
from the effect of the decay products of on BBN [0 . These studies assume monoenergetic 
injection at energy at time r^f,, and their results can be used safely, since the moduli 
emitted by PBHs decay when non-relativistic. One usually considers production of hadrons 
or photons in decay. The constraint due to hadron injection is in principle very significant 
for > 1 TeV, but it is not obvious that can decay hadronically, and moreover it relies on 
assumptions on the cosmic evolution of helium-3 (see, e.g., Ref. 0), which have now been 



proven uncertain (see, e.g., Ref. ||T3 and references therein). Therefore, in the following. 
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only constraints on photon injection are used; the bounds presented will thus be slightly 
conservative. Holtmann et al. obtain in this case: 

< 10-15 for ~ 100 GeV, 

< 10-1^ for ~ 300 GeV, 

< 5 X 10-13 for ~ 1 TeV, and 

< 5 X 10-1° for ~ 3 TeV. 

The error on the upper limit is a factor ~ 4. It results from the uncertainty in the fudge 
factors that enter the T^{m^) relationship when the constraints of Holtmann et al., given in 
the plane m^Y^-r^, are translated into the plane Y^-m^. Note that these constraints assume 
that (p decays into photons with a branching ratio unity, and should be scaled consequently. 
However these limits should also be strengthened by a factor as high as ~ 50 to avoid ^Li 
overproduction, if one assumes that ^Li has not been destroyed in stars in which it has 
been observed Jl^. This constraint is ignored in what follows, as it relies on yet unproven 



assumptions on stellar evolution; this but makes the above constraints more conservative. 
Finally, the observational upper limit on the amount of /x— distortion in the cosmic microwave 
background implies [|T3|, §: Y^ < 10-i3(m<^/100 GeV)!/^ for 20 GeV < < 500 GeV. 



Note that the most stringent constraint on results from the production of the lightest 
of all moduli-like particles in the theory, whose mass would likely be < few x 100 GeV. 
Overall it seems that Y^ < 10"i^ represents a reasonable generic upper limit from BBN. 
Using Eq. (||), this can be rewritten limit on /3: 

and Y^ denotes the upper limit on Y^. This result does not depend on whether PBHs are 



shrouded in a photosphere, as suggested by Heckler |T6[, since moduli are not expected to 
interact with it due to their gravitational interaction cross-section. On the contrary, other 
astrophysical constraints on j3 for Mbh ^ 10^ g are in principle sensitive to the presence of 
a photosphere, as they rely on the direct emission of photons and quarks [^]. 

If i?— parity holds, the constraint on the LSP mass density fli^sp < 1 today implies: 
f3 <2x 10-i^5-i/2^2£(a/3)(xLSp/0.6)-i(MBH/10^GeV)-i/2(mLSp/100GeV)-i. This con- 
straint has been adapted from the study of Ref. ^ and Eq. (||) above. The fraction of 
luminosity carried away by the LSP is xlsp — 0.6, since each spartner produced by a PBH 
will produce at least one LSP in its decay This LSP constraint on (3 is thus less stringent 
than the moduli constraint, provided at least one modulus of the theory has mass < 1 TeV. 

These results have several implications. First of all, the approximation made in Sec. 2, 
namely f2BH <C 1 at all times is justified. In effect, ^bh = P{i/'tiY^'^ at time t in a radiation- 
dominated Universe, since PBHs behave as non-relativistic matter, and therefore at time 
tbr, ^bh — 2.3 X 101''51/^(Mbh/10^ g)/3. Consequently, if P verifies the above upper limits, 
indeed ^Ibk <^ 1 at all times. However, since Eq. is not valid if ^bh = 1 at some time 
t* < tbh, one needs to consider this case as well. 

An order of magnitude of Y^ in this case can be obtained as follows. If ^ tbh, the 
radiation present subsequent to PBH evaporation has been produced in the evaporation 
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process itself. Assuming total conversion of the PBH mass Mbh at time tbh into particles 
(moduli and radiation) of energy aT^n, one finds ^ x^pbr/qTbr, and s ~ (4/3)pbh/^rh, 
where pen is the PBH energy density at evaporation, Trh ~ 2MeV(7io^^'^(^BH/10^g)~^/^ 
is the reheating temperature, and gio = g*/10. Therefore ~ 3 x 10~^°5'io^^^(^</>/6 x 
10-3)(a/3)-l(MBH/10^g)-^/^ well above the previous limits. Note that one would naively 
expect Yfj, ~ since x^jj is the fraction of PBH luminosity carried away by particles. 
However the photons emitted by PBHs carry high energy ~ oTbh and small number density 
~ Pbr/o^Tbr, and their thermalization leads to many soft photons carrying high entropy. 
Nevertheless, this discussion shows that PBHs of any mass should never come to dominate 
the energy density; if this were to happen, PBHs with Mbh ^ 10^ g would produce too 
many moduh, while the evaporation of PBHs with Mbh ^ 10^ g would lead to too low a 
reheating temperature. In particular, scenarios of reheating of the post-inflationary Universe 



by black hole evaporation, as put forward, e.g., in Ref. [0, are forbidden. This result was 
also envisaged in Ref. [Q. 

Finally, the present constraints on (3 exclude the possibility of generating the baryon 
asymmetry of the Universe through PBH evaporation. Indeed Barrow et al. have per- 



formed a detailed computation of the baryon number to entropy density ratio n\^/ s produced 
in PBH evaporation, and find: rib/s ~ 7 x 10^e(a;H/0.01)5'2oo''^(^BH/10^ g)^/^5^/^/5, where e 
is the baryon violation parameter, defined as the net baryon number created in each baryon- 
violating boson decay, xh is the fraction of PBH luminosity carried away by such bosons, 
and other notations are as above. Unless all moduli-like particles are heavier than ~ 3 TeV, 
and i?— parity does not hold, the above constraints on [5 imply rzb/s < 10~^^e, which does 
not suffice since BBN indicates nb/s~4 — 7x 10^^^. 
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